The anomaly in the cosmic-ray positron spectrum 
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Abstract. A recent analysis of cosmic-ray data from a space borne experiment by the AMS col- 
laboration supports the observation of an excess in the cosmic-ray positron spectrum by previous 
balloon experiments. The combination of the various experimental data establishes a deviation 
from the expected background with a significance of more than four standard deviations. The 
observed change in the spectral index cannot be explained without introducing a new source of 
positrons. When interpreted within the MSSM a consistent description of the antiproton spectrum, 
the diffuse gamma-ray flux and the positron fraction is obtained which is compatible with all other 
experimental data, including recent WMAP data. 
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1 Introduction 

Among the cosmic-ray species, antiparticles and dif- 
fuse 7-rays are of particular interest because they are 
produced secondarily in hadronic interactions of pro- 
tons and nuclei with the interstellar medium at low 
rates. Their small abundance makes them a sensitive 
probe for the existence of additional - and possibly 
exotic - cosmic-ray sources which would be visible as 
an excess of particles above conventional expectations. 

One of the most important unsolved questions in 
modern cosmology is the nature of dark matter. The 
most promising dark matter candidate is the weakly 
interacting lightest neutralino, Xii predicted by super- 
symmetric extensions to the standard model of particle 
physics. The annihilation of ncutralinos might consti- 
tute an additional primary source of particles with a 
unique spectral shape which would be determined by 
the parameters of supersymmetry, allowing to put con- 
straints on new physics beyond the standard model. 

A recent reanalysis of the data from the AMS-01 
spectrometer |T] supports the observation of an excess 
of cosmic-ray positrons by the HEAT experiments [5] . 
In this work, we discuss the combined results on the 
cosmic-ray positron fraction e + /(e + + e _ ). Assuming 
that dark matter is largely constituted by neutrali- 
nos, we determine the cosmic-ray preferred parameter 
space of the minimal supersymmetric standard model 
(MSSM) from a simultaneous fit to the cosmic-ray 
positron, antiproton and diffuse 7-ray data. 

2 Cosmic-ray particle propagation 

The public GALPROP code [3] has been used to model 
cosmic-ray particle propagation and calculate the par- 
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ticle spectra as observed near Earth. GALPROP solves 
the propagation equation in a diffusion model with a 
given source distribution for all cosmic-ray species and 
includes convection, diffuse reacceleration, energy loss, 
fragmentation and decay in the interstellar medium. 
The injection spectra of nuclei and electrons before 
propagation are assumed to be power laws in momen- 
tum, and their spectral indices, 7 S and j e , respectively, 
are chosen such that the model reproduces the most 
recent cosmic-ray flux measurements. 

From fitting the propagation model to electron and 
proton flux data we find the most probable values 
of these indices to be 7 S = 2.35 ± 0.03 and 7 e = 
2.50 ± 0.04. In order to determine their errors and 
thus estimate the uncertainties of the model predic- 
tions, the indices have been varied over small intervals 
around their most probable values and the resulting 
predicted fluxes have been compared to the data. The 
\ 2 calculated from the deviation of the data from the 
respective prediction gives the la errors of the injec- 
tion spectral indices. 

Fig. [T] shows the calculated fluxes of electrons and 
protons which are in excellent agreement with the ex- 
perimental data over large energy intervals. The un- 
certainties of the propagation model - the fluxes cal- 
culated with the injection indices at their error limits 
- are denoted by the yellow areas. Below energies of 
several GeV, the individual measurements differ from 
each other due to the time-dependent effect of solar 
modulation. 



3 The cosmic-ray positron fraction 

The challenge of cosmic-ray positron measurements is 
the rejection of the vast proton background. A number 
of balloon borne experiments have delivered positron 
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Fig. 1. a) Compilation of cosmic-ray electron flux 
data from CAPRICE g], HEAT-e* 0, AMS-01 [5] and 
Kobayashi et al. [6]. b) Proton flux data from AMS-01 [7] 
and BESS 8,9 . In both panels, the solid line denotes our 
GALPROP model and its uncertainty. 
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Fig. 2. a) Compilation of recent cosmic-ray positron frac- 
tion data from: AMS-01 (2000) [TT], the AMS-01 reanaly- 
sis p], HEAT-e* and HEAT-pbar 0, CAPRICE g] and 
TS93 10 . b) The combined data together with the back- 
ground model (thick solid line) and its uncertainty (dashed 
lines) as well as the neutralino annihilation signal and sig- 
nal+background for the best fit parameter set (thin solid 
lines). The dotted lines denote the propagation uncertain- 
ties of the signal contribution. 



flux data in the energy range from 0.5 to 50 GeV, 
such as HEAT-e ± and HEAT-pbar 0, CAPRICE gj 
and TS93 [10 . Additionally, the AMS-01 spectrometer 
has measured the positron flux up to 3 GeV |[1 1J in a 
low Earth orbit. In order to extend the sensitivity of 
AMS-01 to energies of up to 50 GeV, a reanalysis of 
the data has been conducted pQ using the conversion 
of bremsstrahlung photons from positrons to achieve 
a proton background suppression of more than 10 5 . 
The result is shown in panel a) of Fig. [2]together with 
previous data. 

In order to simplify data handling, the measure- 
ments on the positron fraction displayed in Fig. [2] a) 
have been combined into one single data set with re- 
gard to asymmetric statistical and systematic errors. 
Details of this procedure as well as a result table are 
given in Ref. [T2]. Panel b) of Fig. |] shows the com- 
bined data together with the model prediction. Above 
energies of 6 GeV, the data exhibit a change in the 
spectral index of positrons which is clearly incompati- 



ble with the expectation for purely secondary positron 
production. Taking into account experimental errors 
as well as the model uncertainty, the significance of 
the deviation amounts to more than four standard de- 
viations. There is no set of propagation parameters 
based on which the GALPROP model would match 
the data satisfactorily. Consequently, the excess in the 
positron flux cannot be explained by the current prop- 
agation models and thus requires a new primary source 
of positrons. 

4 The spectra of cosmic-ray antiprotons 
and diffuse 7-rays 

Using the same procedure as stated in § [3] measure- 
ments of the cosmic-ray antiproton flux from AMS-01 
H31, BESS97 [H|, BESSOO [15] and CAPRICE [16] 
have been combined into one single data set. The re- 
sult is displayed in panel a) of Fig. [3] together with 
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Fig. 3. a) The combined antiproton flux data, b) The dif- 
fuse 7-ray flux data from EGRET [171115] . Both panels: 
modeled background (thick solid line) and its uncertainty 
(dashed lines) as well as the neutralino annihilation sig- 
nal and signal+background for the best fit parameter set 
(thin solid lines). The dotted lines denote the propagation 
uncertainties of the signal contribution. 



the GALPROP calculation. Within the experimental 
errors, the combined data are well in agreement with 
the expectation for purely secondary antiproton pro- 
duction. 

Fig. [3] b) shows the flux of diffuse 7-rays mea- 
sured by the EGRET [171118] experiment. Above en- 
ergies of 1 GeV, the data exhibit a significant excess 
with respect to the model calculation for purely sec- 
ondary production which has been interpreted as an 
additional primary source of 7-rays from neutralino 
annihilations [19 . However, the particular model was 
claimed to be in conflict with the observed antiproton 
spectrum [20]. It has recently been pointed out that 
the excess could also be an artifact from energy mis- 
calibration of the experiment [21] . The discrepancy is 
in principle resolvable by fine-tuning the propagation 
model parameters, which however results in predic- 
tions for other particle species' spectra which are in- 
compatible with experimental data [22] , 



5 Interpretation of the cosmic-ray spectra 
within the MSSM 

5.1 Constraints on the MSSM parameter space 

Measurements of several quantities are used to con- 
strain the parameter space of the MSSM, such as the 
dark matter relic density from WMAP [23 and the 
branching ratios of the rare decays b — > sj [24] or 
Bs — > /i/i [25] , Additional constraints come from the 
LEP2 experiments as lower limits on the neutralio [20] 
and neutral Higgs boson masses |27j . Furthermore, 
measurements of the anomalous magnetic moment of 
the muon [25J suggest low values of the MSSM param- 
eters mo and m 1( / 2 . Fig. [3J shows the plane spanned by 
mo and m 1 / 2 for tan/3 = 40, A = and sign fx = +1 
together with the respective 2a limits derived from the 
above constraints. 



5.2 MSSM parameter scan with cosmic-ray data 

In order to put further constraints on the MSSM pa- 
rameter space from cosmic-ray data, we have conducted 
scans of the plane spanned by the parameters Too and 
TO4/2 for particular fixed values of tan/3. For each of 
the sample points in the plane the contributions to the 
positron fraction and the antiproton and 7-ray spectra 
from neutralino annihilation after galactic propagation 
have been calculated and simultaneously fitted to the 
experimental data together with the GALPROP mod- 
els for the purely secondary background components. 
For the calculations, an isothermal dark matter halo 
profile with a local density of po = 0.3 GeV/cm 3 has 
been assumed. They were performed using the public 
DarkSUSY 4.1 [20], FeynHiggs 1.2.2 [30] and ISAJET 
7.75 [31] packages with the top-quark mass fixed to 
TO tO p=170.9 GeV, A = and signp, = +1. 

The particle fluxes <Pi as observed near Earth can 
be described by adding the calculated signal contri- 
butions, Si, for the individual particle species i to the 
respective GALPROP background models, Bi, accord- 
ing to <Pi = Bi + fi ■ Si. Here, the /, denote individual 
boost factors to allow for a signal enhancement due to 
a possible clumpy nature of the dark matter distribu- 
tion in the solar neighborhood. In this case, we expect 
the individual boost factors to differ significantly from 
each other due to the different travel paths of the parti- 
cle species which are determined by their mean energy 
loss. In particular, the boost factor for the antipro- 
ton signal should be small with respect to the oth- 
ers, since the low synchrotron radiation level of heavy 
particles allows them to be measured almost indepen- 
dently from their production location in the galaxy. 
The boost factors were determined as free parameters 
in fits of the <Pi to the experimental data described in 
§[3] and Sjl 

Fig. Wlshows the combined \ 2 from the simultane- 
ous fits as a function of Too and mi/2 for tan/3 = 40. 
Apparently, the cosmic-ray data clearly favor the fo- 
cus point region at large values of m , and we find 
the best fit parameters to be mo = 1230 GeV and 
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Fig. 4. The plane spanned by the MSSM parameters mo 
and m-t/2 for tan/? = 40, ^0 = and sign/i = +1. Cur- 
rent experimental constraints (2a) are denoted by the solid, 
dotted and dash-dotted lines. The color scale gives the x 2 
of the MSSM fit to the cosmic-ray data. 



m i/2 = 260 GeV. This point is well in agreement with 
all constraints on the MSSM parameter space stated 
in § |5.1| including recent data from WMAP. 

The contributions from neutralino annihilation to 
the individual particle spectra as well as the signal + 
background curves for the best fit parameter set are 
shown in Fig.[2]b) and[3]in comparison with the exper- 
imental data. With the additional primary cosmic-ray 
component from neutralino annihilation, the experi- 
mental data for the positron fraction and the spectra 
of antiprotons and 7-rays can well be reproduced. The 
combined x 2 turns out to be 28 with 33 degrees of free- 
dom, and the boost factors fi are found to be 85 ± 15 
for positrons, 1 ± 0.5 in the case of antiprotons and 
310 ± 50 for 7-rays. 

In the region of the parameter space preferred by 
the cosmic-ray data, neutralinos have a significant hig- 
gsino component of more than 30 % and dominantly 
annihilate into W-boson pairs via t-channel exchange 
of charginos. For the best fit parameters, we find the 
mass of the x? to be 91 GeV and a value for the mass 
of the lightest Higgs boson of 113.7 GeV. 

5.3 Dependence on tan/3 and mtop 

The choice of tan/3 is critical to constrain the MSSM 
parameter space with cosmic-ray data. For varying val- 
ues of tan /3, the combined fits favor a neutralino mass 
between 80 GeV and 120 GeV. Unless tan /3 is higher 
than 50, we always find an overlap of the parameter 
space favored by cosmic rays with the relic density con- 
straints from WMAP in the focus point region. Fur- 
thermore, the preference of cosmic rays in terms of the 
MSSM parameter space is sensitive to the mass of the 
top-quark, whose value is currently known with a pre- 
cision of 1.8 GeV[3"2"]. In particular, for low tan/3 and 



values of m-top > 173 GeV, the focus point region is 
not available unless mo is larger than about 3 TeV. In 
order to put accurate constraints on the MSSM param- 
eter space from cosmic-ray data, the impact of varying 
values of tan (3 and mt op has to be investigated further. 

6 Conclusions 

In this work, the combined recent experimental results 
on the cosmic-ray positron fraction have been pre- 
sented. The data exhibit an excess of positrons above 
energies of 6 GeV which cannot be explained by purely 
secondary positron production alone and thus requires 
an additional primary source of positrons. In this work, 
we interpret this source to be the annihilation of su- 
persymmetric neutralinos constituting dark matter. A 
simultaneous fit to the cosmic-ray positron, antiproton 
and 7-ray data shows that, for particular sets of the 
MSSM parameters, this hypothesis gives a fully con- 
sistent description of the cosmic-ray spectra which is 
compatible with all other experimental data. We find 
that the cosmic-ray data clearly prefer the focus point 
region of the MSSM parameter space but reveal almost 
no sensitivity to tan [3. 
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